The Hidden Landscape of Localization 

Marcel Filoche 1 ' 2 *, Svitlana Mayboroda 3 

O 

Physique de la Matiere Condensee, Ecole Polytechnique, CNRS, 

! ^ 91128 Palaiseau, France 

(N 2 CMLA, ENS Cachan, CNRS, UniverSud, 

i — i 61 Avenue du President Wilson, F-94230 Cachan, France 

3 Department of Mathematics, Purdue University, 

150 North University Street, West Lafayette, Indiana 47907-2067, USA 

*Marcel Filoche; E-mail: marcel.filoche@polytechnique.edu. 

s 

I 

Wave localization occurs in all types of vibrating systems, in acoustics, me- 
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chanics, optics, or quantum physics. It arises either in systems of irregular ge- 



j> ometry (weak localization) or in disordered systems (Anderson localization). 

We present here a general theory that explains how the system geometry and 

O 



the wave operator interplay to give rise to a "landscape" that splits the system 
into weakly coupled subregions, and how these regions shape the spatial dis- 
tribution of the vibrational eigenmodes. This theory holds in any dimension, 
^ for any domain shape, and for all operators deriving from an energy form. It 

b 

encompasses both weak and Anderson localizations in the same mathematical 
frame and shows, in particular, that Anderson localization can be understood 
as a special case of weak localization in a very rough landscape. 



Introduction Every object in nature vibrates. The vibrations can be acoustic waves in a 
medium, mechanical deformations of a rigid plate, electromagnetic waves in a cavity or quan- 
tum states such as the electronic states in a crystal. In all these cases, the linear behavior of 
the system can be reduced to the knowledge of the vibrational modes, i.e., the eigenfunctions 
and eigenvalues of the spatial differential operator associated to the wave equation. A stun- 
ning property exhibited by these eigenfunctions in irregular or disordered systems is known 
as localization: in some cases, the mode amplitude is small everywhere except in a very lim- 
ited subregion of the domain, even though there is no clearly visible obstacle preventing the 
vibration to propagate in the rest of the domain, away from its main existence subregion. 

It has been acknowledged that there are several types of localization, each exhibiting a spe- 
cific behavior in space or in frequency. First, when caused by the irregular or complex geometry 
of the vibrating domain, localization is classified as weak. It is characterized by a slow decay 
of the mode amplitude away from its main existence subregion (much slower than exponen- 
tial) (1-4). Secondly, localization can arise due to a disorder quenched in the system (5, 6). This 
phenomenon, called Anderson localization, has fascinated scientists since its discovery in 1958 
and spurred a wealth of literature (7). In that case, the mode amplitude decays exponentially 
away from its main existence subregion and the localization is classified as strong. Finally, a 
third and quite different type of localization occurs at high frequencies in specific domains that 
possess stable orbits such as whispering gallery or bouncing ball modes (8). Confined modes 
known as scar modes then appear asymptotically localized near these orbits. 

Until today, there has been no theory able to explain how the geometry of the domain or the 
nature of the disorder is related to the localization of vibrations, to predict in which subregions 
one can expect localized eigenmodes to appear, and in which frequency range. Moreover, the 
question of whether the different types of localization are linked remains open. Consider, for 
instance, the system depicted in Fig. [T] (left). It has a non trivial shape, possesses 2 inner blocked 



points (in the left upper region), one crack on the right upper boundary and a bottleneck in its 
lower region. This can represent either a flexible membrane of complex shape (the spatial differ- 
ential operator is then the Laplacian), or a rigid thin plate (the operator being the bi-Laplacian). 
From the knowledge of this geometry only, how can one determine whether and where to expect 
localization in this structure? In the present paper, we show that one theory can answer these 
questions. It unifies weak and Anderson localization, and reveals inside each system a hidden 
landscape that determines the localization subregions, the strength of the confinement, and its 
frequency dependency. 

The landscape mapping In mathematical terms, a vibrating system is governed by the wave 
equation associated to a suitable elliptic differential operator L. The latter is determined by 
the nature of vibration and the medium. For instance, the Laplacian L = — A describes the 
vibration of a soft membrane in 2D, acoustic waves, the quantum states inside a box or cavity; 
variable coefficient second order operators L = — div (A(x) V) pertain to the aforementioned 
phenomena in inhomogeneous media; and the bi-Laplacian A 2 addresses thin plate vibrations 
in 2D. A study of the vibrational properties of the system can be reduced to the investigation of 
the eigenmodes of L defined by 

Lip = Ay? in Q, (p\ dn = (1) 

We demonstrate here that all eigenmodes are controlled by the same function which has a 
decisive impact on localization properties. To be specific, for every A and ip as above 

\<p(x)\<\ [ \G(x,y)\dy, ViGfi, (2) 

Jn 

where the mode y? is normalized so that sup Q \ip\ = 1, and G(x,y) is the Green's function 

solving L G(x, y) = 5g(y) with zero data on the boundary. (For the sake of brevity, we assume 

Dirichlet boundary data. Other types of boundary conditions will be addressed in forthcoming 
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publications). In other words, there exists a function u(x) = J n \G(x, y)\ dy, independent of 
the eigenmode ip such that 

\<p(x)\ < \u(x) (3) 

Recall that the Green's function is positive for the Laplacian and, more generally, for all second 
order differential operators. In that case u admits a remarkably simple definition. It is the 
solution to the Dirichlet problem 

L« = 1 in H, u\ga = (4) 

For a higher order operator the boundary condition u\qq = in ([T]) and Q refers to the natural 
Dirichlet data, including vanishing derivatives of appropriate orders. We defer rigorous math- 
ematical justification of ©-Q to (9). In physical terms, u can be interpreted, for instance, as 
the steady-state deformation of a membrane under a uniform load. 

Through inequality ([3]), the landscape u compels the eigenmodes to be small along its lines 
of local minima (called valleys throughout the paper). As we will show, the network of these 
valleys, a priori invisible when looking at the domain, but clearly identifiable on the graph of 
u, operates as a driving force that determines the confinement properties for both weak and 
Anderson localization. 

Let us illustrate this with an example. Figure[T]displays two maps of u computed in the same 
complex geometry for the Laplacian (center) and the bi-Laplacian (right), respectively. The 
streamlines (lines of the gradient) have been plotted in thin black in order to clearly pinpoint the 
valleys, which are highlighted as thick red lines. The two cases expose dramatically different 
patterns. For the Laplacian, one can observe two valleys, and only one of them splits the domain 
into two disjoint subregions. In the bi-Laplacian case, five valleys form a network yielding a 
partition of the domain into four disjoint subregions. 

We now show rigorously that for any domain and any differential operator the network of 




—0.025 
n no 



x10 



0.015 



0.01 



0.005 




Figure 1: Left: Geometry of a complex domain, with a bottleneck (in the lower part), two 
inner blocked points, and an inward crack (on the right upper side). The question is: is there 
localization in this structure, which modes will be localized and where? Center and right: 2D 
representations of the landscape u for the Laplacian (center) and the bi-Laplacian (right) in the 
same domain. The colors correspond to the height of the landscape. The (thin black) streamlines 
help to detect the valley lines, which are then highlighted as thick red curves. These valleys 
delimit 2 subregions of localization for the Laplacian and 4 subregions for the bi-Laplacian. 

valleys built as above triggers wave localization by effectively separating the original domain 
into weakly coupled vibrating regions. 

The formation of localized modes Consider a subregion Vt\ of the original domain Vt carved 
out by the valley lines. It can be thought of, for instance, as any of the four subregions in 
Figure [Tj right. By construction, u is relatively small (locally minimal) along the boundary 
of f^i. Thus, for low eigenvalues A, inequality ([3]) provides a severe constraint on the mode 
amplitude y? along the boundary dVt\. As a result, any eigenmode ip of the entire domain can 
locally be viewed as a solution to the problem 

L ip = X ip in Qi, (5) 

ip = on dQ.1 fl dSl, and cp = e on d£l\ \ dil, (6) 



where e(x) is a quantity smaller than Xu(x) on the boundary of fii. Observe that boundary 
value problem (J?])-© is, in fact, akin to the eigenvalue problem in the subregion Qi alone: 



the differential equation inside the subregion is identical, but on the boundary ip is small in (|6]) 
rather than just being zero as an eigenvalue problem on VL\ would normally warrant. 

Using properties of the resolvent (L — A/) -1 and spectral decomposition, one can establish 
the following estimate which plays an essential role in understanding the origin of localization: 

»*<*> £ t 1 + dim) w <7) 

See (9) for the proof. Here, dn^X) is the distance from A to the spectrum of the operator L in 
the subregion f^ (defined as: dn 1 (X) = min {|A — A^qJ}, the minimum being taken over all 
eigenvalues (A^nJ of L in fix), and ||e|| is the L 2 -norm of the solution to Lv = in fii with 
data e on dVti (in the sense of (|6])). In particular, ||e|| becomes arbitrarily small as e in Eq. (|6]) 
vanishes. 

The presence of da,i (A) in the denominator of the right-hand side of Eq. |7]) assures that 
whenever A is far from any eigenvalue of L in fix in relative value, the norm of cp in the entire 
subregion, ||<^||z,2( ni ), has to be smaller than 2||e||. Consequently, such a mode cp is expelled 
from Qi and must "live" in its complement, exhibiting weak localization. Conversely, tp can 
only be substantial in the subregion Vti when A almost coincides with one of local eigenvalues 
of the operator L in f^. Moreover, in that case Eq. ([5]) yields the conclusion that <p itself almost 
coincides with the corresponding eigenmode of the subregion f^i. 

Thus, we obtain a rigorous scheme elucidating the formation of weak localization. In any 
subregion delimited by the valleys of u, an eigenmode of VL has only two possible choices: (1) 
either its amplitude is very small throughout this subregion, or (2) this mode mimics (both in 
frequency and in shape) one of the subregion's own eigenmodes. Consequently, a low frequency 
eigenmode can cross the boundary between two adjacent subregions only if they possess two 
similar local eigenvalues. More generally, & fully delocalized eigenmode can only emerge as a 
collection of local eigenmodes of all subdomains when they all share a common eigenvalue. 
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Figure 2: two localized modes of the Laplace operator (left) and four localized modes of the 
bi-Laplace operator (right). In both cases, the valley network (displayed in red) obtained in 
Figure [TJ accurately predicts the number and the location of the localization subregions. 



The role played by the valleys of u in localization is even more astonishing when exploring 
the examples. Let us revisit the domain in Figure[TJ Figure[2j left, displays two localized modes 
(number 1 and 5) of the Laplacian plotted together with the valley lines computed in Figure [TJ 
center. Not only these modes obey the pattern predicted by the landscape u, but there exists no 
eigenmode confined to a smaller subregion. In Figure [2j right, modes 1,2, 4, and 6 of the bi- 
Laplacian are displayed in the same domain together with the valley lines from Figure [TJ right. 
We observe here a very different and much larger variety of localization behaviors. However, 
once again, the location, the shape, and the number of the localization subregions exactly match 
the partition of the domain generated by the valley network of the landscape of u computed in 
Figure [TJ 



The effective valley network and Weyl's law Not only the location of the subregions is 
revealed by the network of valleys, but the height and shape of the valleys gives access to statis- 
tical information on the evolution of the localization properties in all frequency ranges. While 
geometrically the network of valleys is determined by u and does not depend on a particular 
eigenmode, the strength of the confinement of an eigenmode dictated by Eq. ([3]) diminishes as 
A grows. Indeed, given the normalization chosen in Eq. (|2]), Eq. ([3]) represents an effective con- 
straint only at those points {x} where A u(x) is smaller than 1. In other words, an eigenmode 
can actually see only the portion of the initial landscape which satisfies u < 1/A. We refer to 
the latter as the effective valley network. 

By definition, at relatively low frequencies the effective network is identical to the full one. 
However, as A increases, the effective network progressively disappears. Subregions that were 
initially disjoint begin to merge to form larger subregions. Above a critical value of A, what is 
left of the effective network allows one subregion to percolate throughout the system: starting 
from that value of A, completely new fully delocalized modes can appear. 

Note that the precise quantitative information on evolution of the effective network with 
the growth of A is already encoded in the original mapping of u. In this vein, it would be 
extremely useful to estimate the underlying rate of growth of the eigenvalues. Unfortunately, 
in the full generality of our set-up the Weyl's law is not available yet. Roughly speaking, it 
predicts that for any given A > the number of the eigenvalues of L below A is asymptotically 
\ d i 2m 5 where d stands for the dimension and 2m is the order of the differential operator L. This 
estimate is widely used in practice and can be employed in the present context to determine 
the range of frequencies with high response to the impact of u, that is, the range of thoroughly 
localized eigenmodes. In this context, one has to stress that in a domain of self-similar boundary, 
smaller copies of a valley line would appear at all scales, triggering mode localization for an 
infinite number of eigenvalues. In that case, even though weak localization essentially affects 
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low frequency eigenmodes, one would always find localized modes in the high frequency limit 
governed by Weyl's law. 

From weak to strong: Anderson localization Since Anderson's work in 1958 (5), strong lo- 
calization is known to arise due to the presence of structural disorder in a system. For instance, 
in a gas of non interacting electrons in a crystal, this phenomenon can induce a transition be- 
tween metallic and insulating behavior as the electronic states become highly localized (6, 10). 
More recently, disorder-induced localization has been demonstrated to be a very general phe- 
nomenon also observed in acoustics {11), in microwaves (12), or in optics (13, 14). 

We present here a fundamentally novel approach to Anderson localization, resting upon 
the theory developed in the previous sections. To this end, we compute eigenmodes of the 
Schrodinger operator with a random potential modeled as follows. The original domain is cho- 
sen to be a simple unit square. It is divided into 400 = 20 x 20 smaller squares. On each of these 
smaller squares, the potential V is constant, its value being determined at random uniformly be- 
tween and Vmax (here V m3iX = 8000, see Figure [3j top left, the arbitrary energy units are taken 
such that h 2 /2m = 1). The corresponding Schrodinger operator H = —A + V is therefore a 
second order elliptic operator with variable coefficients falling under the scope of our theory. 
The valley landscape is further obtained by solving H u = 1 in the entire square, with Dirichlet 
boundary conditions on 4 sides (see Figure [3} top right). One can observe a complex relief, 
characterized by a network of interconnected valleys of varying depths. Using the streamlines 
as guidelines, the valley lines are drawn in Figure [3j bottom. The thick lines represent the deep- 
est valleys, while the moderately thick lines correspond to the shallower ones. This intricate 
network reveals a complex partition of the domain into a large number of subregions that was 
impossible to guess by just looking at the random potential at hand (cf. Figure [3j top left). In 
the continuous limit of a Gaussian random potential, one can conjecture that this network would 
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Figure 3: Top left: random potential V(x). The 2D square domain is divided in 20 x 20 small 
squares. In each square, the potential V is assigned an independent random value uniformly 
distributed between and V max (here 8000). Top right: 3D view of the landscape of u, obtained 
by solving Lu = 1, where L = —A + V(x). Bottom: 2D color representation of the map of 
u, together with the streamlines. The thicker lines correspond to the deepest valleys and the 
moderately thick lines to the more shallow ones. This landscape draws an intricate network 
of interconnected valleys. One can conjecture that in the limit of a Brownian potential, this 
network becomes scale invariant and exhibit fractal properties. 
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show similar patterns at all possible valley depths, and hence, exhibit fractal properties. 

Figure [4] displays 2D color representations of the amplitude for a number of modes (or 
quantum states), at lower and higher frequencies. For each mode, the effective network of valley 
lines is plotted on top of the amplitude. Since in the present context the eigenvalue is equal to 
the energy, the effective network can be viewed as a subset of the initial landscape of valleys 
subject to constraint u < l/E, E being the mode energy. It is striking to observe how all modes 
are clearly shaped by the valley lines. The fundamental and the first excited states (modes 1 
to 8) are localized completely to one of the subregions defined by the network of valley lines. 
At higher energies (mode 45, 48, 70), the effective valley network starts to shrink, opening 
breaches in the shallowest valley lines. One can see that these modes are still localized, but now 
exactly in the much larger subregions defined by the remaining effective network. However, 
there still exist some small subregions in which one can find localized modes weakly coupled to 
the rest of the domain (modes 47 and 71). At even higher energies (modes 97-99), the effective 
valley network is mostly disconnected, allowing a subregion to percolate throughout the entire 
domain: delocalized states appear. 

For sake of simplicity, simulations have been carried out in 2D. In 3D, the valleys are not 
lines but surfaces of minimal value of u. The entire valley network has the shape of a foam 
which separates the domain into a large number of subregions. When the energy increases, the 
effective valley network evolves by opening gaps in the walls separating adjacent subregions. 
The localization length then increases accordingly as the average size of a subregion. 

One can examine in detail the strength of the localization by plotting the mode amplitude on 
a logarithmic scale (Figure [5]). By doing so, one can notice that the level curves are on average 
equally spaced, which corresponds to an exponential decay away from the existence subregion. 
Even more precisely, we once again observe our landscape at work: the amplitude of the mode 
appears to stay within the same order of magnitude inside each subregion, the decay essentially 
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Figure 4: Spatial distribution of several quantum states in the random potential of Fig. [3} On 
top of each state is drawn the effective valley network corresponding to the state energy. One 
can clearly observe that all modes are localized exactly in one of the subregions delimited by 
the effective valley network. 
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Figure 5: Logarithmic plot (in log 10) of the amplitude for two excited states in the potential 
shown in Fig. [3j One can observe here in detail how strong localization emerges from weak 
localization. Both modes are located in one of the subregions delimited by the valleys ways. 
Moreover, their decay is shaped by neighboring subregions. Firstly, the mode amplitude is more 
or less uniform within any subregion. Secondly, when going away from the main subregion, 
crossing a valley line corresponds to a decrease by one or two orders of magnitude. This leads 
to an exponential decay for distance larger than the typical size a subregion. 

occurs when crossing a boundary between two adjacent subregions (this boundary corresponds 
to a valley line). This is particularly clear for both modes in Figure [5] in which their principal 
existence subregion appears dark red, all nearest neighboring subregions appear light red, the 
second nearest neighboring subregions essentially orange, etc. Therefore, each subregion is 
weakly coupled to its neighbors, the mode decaying by a somewhat constant factor each time 
it crosses a valley line away from the center subregion of existence. Hence, when zoomed 
towards any subregion, mode localization is of the weak type. However, due to the intricacy of 
the valley network, the succession of regularly spaced valley lines in the effective network yields 
an exponential decay of the mode away from its center subregion, over distances much larger 
than the typical size of a subregion. As a consequence, the mode appears to be exponentially 
confined: strong localization emerges from successive weak localizations. 

Therefore, the exact same concepts and the exact same theory both account for localization 
in a domain of irregular geometry and localization induced by the presence of disorder. More- 
over the notion of effective network extremely accurately predicts the shape of the localization 
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subregions at any energy and therefore allows us to understand the apparition of delocalized 
states above a critical energy: this critical value is the smallest energy for which the effective 
network opens all valley lines and allows one subregion to percolate throughout the system. 
Remarkably, computing the landscape of u and the effective valley network at any energy does 
not require any knowledge on the quantum states of the system, but yet gives access to accu- 
rate information on the confinement of the states within any energy range. One also should 
add that, although the landscape u is obtained by simply solving one linear system, it neverthe- 
less depends in a complex way upon the quenched disorder introduced in the potential V. In 
particular, different types of randomness would lead to different localization properties. More 
generally, the question of localization in any disordered system can now be mathematically re- 
formulated in the simple following way: How does the valley network of the landscape of u (its 
connectivity, sizes, and height), depend on the properties of the disorder in VI 

Localization also directly affects the transport properties in a random medium (15): at zero 
temperature, the system is insulating if all the quantum states occupied by the electrons are 
localized. This condition can be reformulated in terms of the landscape: for the system to be 
insulating at T=0, its effective valley network computed at the Fermi level (in other words, the 
network of valley lines for which u < 1/Ep) must allow one subregion to percolate through the 
entire system. In 2D, if a subregion crosses the entire domain then its effective valley network 
has to be disconnected. The percolation transition of the subregion thus exactly corresponds to 
the percolation transition of the effective valley network itself. 

This theory also allows us to understand how temperature can modify the conductance of 
the system. On one hand, its increases the amplitude of disorder that enters the Schrodinger 
equation. This tends to localize the quantum states and to reduce the conductance. On the 
other hand, it explores excited states at higher energies, which are less localized, and it couples 
the states, allowing electrons to hop between states. The existence of a transition at finite 
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temperature in a macroscopic system therefore depends on the outcome of the competition 

between these effects. 

Finally one should recall that this theory is formally valid in any dimension. In particular, 

it means that the localization or derealization properties of a system of N correlated electrons 

can also be analyzed, in theory, through the knowledge of the statistical properties of the valley 

/ n \-i 

network for the landscape u — I — 2_. ^i + V J 1 in a 3iV-diniensional space. 
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Conclusion What emerges here is a new picture of localization. There are, in fact, not three 
but only two types of localization: low frequency localization, described by the landscape theory 
developed in this paper, and scar or high frequency localization, predicted by the stable orbits 
of the domain. 

Our findings demonstrate that low frequency localization is a universal phenomenon, ob- 
served for any type of vibration governed by a spatial differential operator L that derives from 
an energy form. The geometry of the domain and the properties of the operator interplay to 
create a landscape u which entirely determines the localization properties of the system. First, 
its network of valley lines (in 2D) or surfaces (in 3D) creates a partition of the initial domain 
into disjoint subregions which shape the spatial distributions of the vibrational modes and iden- 
tify precisely the regions confining vibrations. Second, the depth of these valleys determines 
the strength of the confinement within each subregion. The localization of a given mode of 
eigenvalue A is controlled by the effective valley network defined as the subset of the entire 
original valley network subject to the condition u < 1/A. Consequently, the relative number 
of localized modes decreases at high frequencies. After partitioning by the landscape, a com- 
plex vibrational system can be understood as a collection of weakly coupled subregions whose 
coupling increases with frequency. 

The theory holds for systems of irregular geometry as well as for disordered ones. In this 
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framework, Anderson localization arises as a specific form of weak localization, "strengthened" 
by the extremely rough landscape generated by a random potential. More generally, the macro- 
scopic properties materials or systems in which localization plays an essential role, can now be 
reformulated from the geometrical and analytical characteristics of the effective valley network. 

This theory of localization opens a number of new problems. Let us mention a few: In the 
case of domain with fractal boundary, can one relate the asymptotic distribution of eigenmodes 
with the scaling properties of the valley network? Can one deduce the thermodynamical be- 
havior of non interacting bosons or fermions in a disordered system from the knowledge of its 
valley network? What are the statistical properties of the landscape of N interacting particles? 

Finally, one should stress that the effective valley network is promising to become a new tool 
of primary importance for designing systems with specific vibrational properties. To this end, 
future studies should investigate in detail the relationship between the geometry of a system 
(irregular or fractal), the characteristics of the wave operator (order, non-homogeneity, possibly 
stochastic), and the properties of one key mapping, the resulting valley network. 
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00-18. The second author is partially supported by the Alfred P. Sloan Fellowship, the NSF 
CAREER Award DMS 1056004, and the US NSF Grant DMS 0758500. The authors also 
wish to thank Guy David for fruitful discussions. 
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SUPPORTING ONLINE MATERIAL 



Materials and methods 



Preliminaries Let f2 be an arbitrary bounded open set in MJ 1 and let L be any elliptic differen- 
tial operator associated to a symmetric positive bilinear form B (an energy integral). Essentially 
all elliptic operators governing wave propagation, whether in acoustics, mechanics, electromag- 
netism, or quantum physics, are associated with an energy and fall into this category. The most 
prominent examples include: 

the Laplacian L = — A, B[u,v] = / VuVvdx, (1) 

Jn 

the Hamiltonian L = -A + V(x), < V(x) < C, B[u, v] = I WuWv + Vuv dl&) 

Jn 

the bilaplacian L = A 2 = -A(-A), B[u,v]= AuAvdx, (3) 

Jn 

and finally, any second order divergence form elliptic operator 

L = -divA(x)V, B[u,v]= A(x)VuVvdx, (4) 

Jn 

where A is an elliptic real symmetric n x n matrix with bounded measurable coefficients, that 

is, 

n 

A(x) = {aa(x)}^ =1 , xeQ, aij e L°°(Q), Yl "ijWtei > MtW ve e R n , (5) 

for some A > 0, and a^ = a^ Vi, j = 1, ..., n. 

In general, L is a differential operator of order 2m, m E N, defined in the weak sense: 

/ Luvdx :— B[u, v], forevery u,v G H m (Q), (6) 

Jn 
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where B is a bounded positive bilinear form and H m (tt) is the Sobolev space of functions given 
by the completion of C£°(Cl) in the norm 

IMI.ff m (n) := \\^ mu \\L 2 (n), (7) 

and V m u denotes the vector of all partial derivatives of u of order m. Recall that the Lax- 
Milgram Lemma ascertains that for every / G (if m (f2))* =: H~ m (Q) the boundary value 
problem 

Lu = f, «6 H m (n), (8) 

has a unique solution understood in the weak sense. 

Remark. The weak solution formalism is necessary to treat the Dirichlet boundary problem on 
an arbitrary bounded domain. When the boundary, the data, and the coefficients of the equation 
are sufficiently smooth, the weak solution coincides with the classical solution and ([8]) can be 
written for a second order operator as 

Lu = f, u\ dn = 0, (9) 

where u\qci is the usual pointwise limit at the boundary. For a 2m-th order operator the deriva- 
tives up to the order m — 1 must vanish as well, e.g., 

A 2 u = f, u\ dn = 0, d v u\ 9 a = 0, (10) 

where d u stays for normal derivative at the boundary. In any context, the condition u e H m (VL) 
automatically prescribes zero Dirichlet boundary data. On rough domains the definitions akin 
to (|9]), ( |T0| ) might not make sense, i.e., a pointwise boundary limit might not exist (the solution 
might be discontinuous at the boundary), and then the Dirichlet data can only be interpreted in 
the sense of©. 

The classification of domains in which all solutions to the Laplace's equation are continuous 
up to the boundary is available due to the celebrated 1924 Wiener criterion (16). Over the years, 
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Wiener test has been extended to a variety of operators. We shall not concentrate on this issue, 
let us just mention the results covering all divergence form second order elliptic operators (IT), 
and the bilaplacian in dimension three (18). 

Here we shall impose no additional restriction on <9f2 or on the coefficients and work in the 
general context of weak solutions. 

For later reference, we also define the Green function of L, as conventionally, by 

L x G(x,y) = 5 y (x), foTallx,yeQ, G(-,y) G H m (n) forallyGfi, (11) 

in the sense of @, so that 

/ L x G{x,y)v(x)dx = v{y), y G £1, (12) 

JM. n 

for every v G H m (Q). It is not difficult to show that for a self-adjoint elliptic operator the Green 
function is symmetric, i.e., G(x, y) = G(y, x), x,y G O. 

Control of the eigenfunctions by the solution to the Dirichlet problem: the landscape Let 

us now turn to the discussion of the eigenfunctions of L. Unless otherwise stated, we assume 
that L is an elliptic operator in the weak sense described above and that the underlying bilinear 
form is symmetric, i.e., that L is self-adjoint. 

The Fredholm theory provides a framework to consider the eigenvalue problem: 

Lip = \(p, <pe # m (ft), (13) 



where A G R. If there exists a non-trivial solution to (13), interpreted, as before, in the weak 



sense then the corresponding A G R is called an eigenvalue and (p G H m (Vi) is an eigenvector. 

Proposition 0.1. Let Q, be an arbitrary bounded open set, Lbe a self-adjoint elliptic operator on 
Q, and assume that if G H m (Q) is an eigenfunction ofL and A is the corresponding eigenvalue, 
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i.e., (p~3|) is satisfied. Then 



\<p(x)\ 



|y||L°°(n) 



< Xu(x), for all x E £1, 



(14) 



provided that ip e L°°(f2), wi(/i 



■u(x) = / |G(x, y)|dy, iGU. 
Jn 



(15) 



7/J m addition, the Green Junction is non-negative (in the sense of distributions), then u is 
the solution of the boundary problem 



Lu = l, ue H m {fl). 



(16) 



Remark. The Green function is positive and eigenfunctions are bounded for the Laplacian ([T]), 
the Hamiltonian (|2]), all second order elliptic operators ([4]) in all dimensions due to the strong 



maximum principle (see, e.g., (79), Section 8.7). Hence, for all such operators ( |T4| ), ( fT6[ ) are 
valid. 

The situation for the higher order PDEs is more subtle. In fact, even for the bilaplacian the 



positivity in general fails, and then one has to operate directly with ( 15 ) 



Proof. By ( |13| ) and ( |TT| ) (with the roles of x and y interchanged) and self-adjointness of L, for 
every x 6 fi 



cp(x) = / (p(y)L y G(x,y)dy= / L y ip(y)G(x,y) dy = / \ip(y) G(x,y)dy, (17) 
Jn Jn Ju 

and hence, 



<p(x)\ < A |MU«( n ) / |G(a;,j/)|dj/, i6fi, 



(18) 



as desired. Moreover, if the Green function is positive, 



/ \G(x,y)\dy= / G(x,y)-ldy, xeQ, 
Jn Jn 



(19) 
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which is by definition a solution of (|T6|). □ 



Vaguely speaking, the inequality ( fl4| ) provides the "landscape of localization", as the map 
of u in ( p"5| ) - p6] ) draws lines separating subdomains which will "host" localized eigenmodes. 
We now discuss in details the situation on these subdomains. 

Analysis of localized modes on the subdomains The gist of the forthcoming discussion is 
that, roughly speaking, a mode of Vt localized to a subdomain D cVt must be fairly close to an 
eigenmode of this subdomain, and an eigenvalue of Vt for which localization takes place, must 
be close to some eigenvalue of D. 

To this end, let f be one of the eigenmodes of fi, which exhibits localization to D, a subdo- 
main of fl This means, in particular, that the boundary values of f on dD are small. In fact, 
in the present context the correct way to interpret "smallness" of if on the boundary of D is in 
terms of the smallness of an L-harmonic function, with the same data as f on dD. 

To be rigorous, let us define s = s^ > as 

e = |H|i,2( D ), where v G H m (D) is such that 

w := ip — v G H m (D) (that is, the boundary values of f and v on dD coincide), (20) 

and Lv = on D in the sense of distributions. 

Then the following Proposition holds. 

Proposition 0.2. Let Q be an arbitrary bounded open set, L be a self-adjoint elliptic operator 
on fl, and f G H m (Q) be an eigenmode of L. Suppose further that D is a subset ofVL and 
denote by e the norm of the boundary data off on dD in the sense of ( [20] ). 

Denote by A the eigenvalue corresponding to if. Then either A is an eigenvalue ofL in D or 



MU 2 (D) < 1+ max 

A fc (£>) 



x _ A*p) 



A 



(21) 
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where the maximum is taken over all eigenvalues of L in D. 



Proof. First of all, note that ( 20 ) implies 



(L — X)w = Xv on D, 



(22) 



as usually, in the sense of distributions. If A is an eigenvalue of D, there is nothing to prove. If 
A is not an eigenvalue of D, we claim that 

1 



MIl 2 (d) — max 



Xv 



L 2 (D)- 



(23) 



A h (fl) [\X-X k (D)\ 

Indeed, in our setup, the eigenvalues of L are real, positive, at most countable, and moreover, 
there exists an orthonormal basis of L 2 (D) formed by the eigenfunctions of L on D, {ipk,D}k- 
In particular, for every / G H m (D) C L 2 (D) we can write 



/ = ^cMA, 



c k (f) 



f ipk dx, 



(24) 



D 



\2W2 



with the convergence in L 2 (D), and ||/||l 2 (d) = (X)& c k{f) 2 ) ±/ *- Moreover, such a series 
Sfc c fe(/)^fc converges in H m (D) as well and {^fc.z?}fc form an orthogonal basis of H m (D). 
These considerations follow from ellipticity and self-adjointness of L in a standard way using 
the machinery of functional analysis (see, e.g., (20), pp. 355-358 treating the case of the second 
order operator of the type (Q). 

Therefore, for every A not belonging to the spectrum of L on D and w G H m (D) with 
(L - X)w G L 2 (D) (cf. ©) we have 



\\(L- X)w\\ L 2 {D) 



y^c fc ((L- X)w)ij)i 



(25) 



L 2 (D) 



where 



c k ((L - X)w) 



(L — X)w^pkdx= / w (L — X)ijjkdx 



D 



D 



(Aa(D)-A) / w^ k dx = (X k (D)-X)c k (w). 
Jd 



(26) 
(27) 
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Hence, 



[L- X)w\\ L 2 



(D) 



)X\k(D) ~ X)c k (w)i) k 



> min \X k (D) - X\ y^c k { 

A fc (D) \ ^ 



1/2 



£(A*(£>) - X) 2 c k ( 

L 2 (D) \ k 

1/2 



W 



W) 



min \X k (D) - X\ ||w||z, 2 (r>), 



which leads to inequality ( |23| ). 
Going further, ([23]) yields 



W £2mi < max 
A fc (D) 



1 



A fc (L>) 



A 



-i 



I ^ II l 2 (di < m & x 



A fc (L>) 



-1 



e, (28) 



and therefore, 



as desired. 



MU 2 (D) < 1 + max 

A fc (Z?) 



A fc (£>) 



(29) 
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